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Abstract 

We show that an R-parity violating supersymmetric scenario which can account for the atmo- 
spheric v^ deficit has testable signals at the Fermilab Tevatron with upgraded energy and luminosity. 
The explanation of neutrino masses and maximal v^~Vt oscillation in terms of bilinear R-violating 
terms in the superpotential associates comparable numbers of muons and tau's resulting from de- 
cays of the lightest neutralino. We show that this should lead to like-sign dimuons and ditaus with 
substantial rates, in a form that separates them from standard model backgrounds and other sig- 
nals of supersymmetry. One here also has the possibility of completely reconstructing the lightest 
neutralino. 

Existing data [||] strongly point towards nondegenerate, massive neutrinos and large mixing between 
v^ and Vr (or, perhaps, a sterile state to which v^ oscillates). To decide on which extension of the 
standard model fits this picture best, one also needs to think of other testable consequences of any 
suggested scenario. In this spirit, here we predict some signals of one such candidate theory 0: a 
supersymmetric (SUSY) scenario [^ where R-parity breaking y] bilinear terms |5|, g mix neutrinos 
and neutralinos Q, and thence lead to massive neutrinos. 

Admitting of lepton number violation, the superpotential for the minimal SUSY standard model 
(MSSM) can be extended to include the following terms: 

W^ = €iLiH2 + XijkLiLjE^ -\- X^jf^LiQjDl (1) 

Here we take the most econonomic aproach and consider only the bilinear terms in the second and 
third generations, viz. €2,3^2,3^2- These terms can be responsible for generating one tree- level neutrino 
mass; also, the retention of the bilinears in the second and third generations ensures f^ — Vr mixing of 
the kind suggested by the atmospheric neutrino data [0] . 



Performing a rotation 



we remove them from the superpotential, and reabsorb them in fj,'H[H2 



(where fi' = {iJ? + e^ + £3)^'^, and H^ is a linear combination of Hi, L2 and L3). But the effects of €2 
and €3 are still manifest; apart from newly produced trilinear terms, they show up through the scalar 
potential, where non-zero sneutrino ueu's u^ and Vj- are in general induced. These vev^s give rise to 
neutrino-gaugino cross terms, that change the structure of the neutralino mass matrix 
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M and M' are the SU(2) and U(l) gaugino masses, and v'l = {H[), V2 = (-^2)- Thence, the state 

z/3 = Uj- cos 9 + v^ sin 9 where tan 9 = U/i/ur (3) 

acquires a tree- level mass. For mz <^ ^', this mass is 



2 M MM' - m\ MJiJ sin 2/3 



(4) 



where 5f^M = g^M' + g' M, and tan /3 = v^Iv'y- The parameter vq = (v^ + u^)"^'^ is a basis-independent 
measure of R-parity violation; if vq ~ 100 keV (say, for //' = —500 GeV, tan/3 = 5 and m^o ~ 100 
GeV) the atmospheric neutrino problem can be explained in terms of oscillations^. The decay of the 
lightest neutralino Xi (assumed to be the LSP-lightest SUSY particle) could provide the crucial test of 
this model: In fact, if enough massive, Xi has two-body decays of the type Xi — ^ ^^ and Xi — ^ ^iZ; 
1 = H,T (Roy-Mukhopadhyaya in g ; ||, |, |ll| ) . 

The production and cascade decays of superparticles at a hadronic collider (till the LSP is reached) 
is expected to be controlled by R-parity conserving interactions. Thus all SUSY processes at the 
Tevatron will end up in Xi pairs. Two-body decays of a neutralino will lead often to lepton-T^ final 
states; a conservative estimate predicts B(xi -^ IJ'W) ~ B(xi — > tW) ~ 35% for maximal mixing angle 
[^. Half of these events contain like-sign dileptons (LSD) due to the Majorana character of neutralinos 
[^]. It is this real W + LSD signal (with correlated numbers of dimuons and ditaus, but a depletion 
of dielectrons) that can help us testing the proposed model. Such a signal is of particular interest for 
the CDF detector, sensitive to LSD's; it has already been used to look for R-parity violation of other 
types with 107 pb^^ data sample obtained at ^/s = 1.8 TeV [^]. Here we present our results for the 
upgraded Tevatron running at ^/s = 2 TeV, with a luminosity of 10 fb~^. 

We assumed that the only QCD process that leads to LSP pairs is the pair-production of squarks 



(of five degenerate flavours) whose lower mass limit is approaching 300 GeV [14|. The effective de- 
coupling of the gluino can be indeed derived from the condition that the lightest neutralino heavier 
than the W, assuming gaugino mass unification (that we use to reduce the number of parameters). 
We retained all possible cascade channels in our calculation, and took particular note of the small but 
non-negligible contributions from X2X1 and xtxi {X2 ^^^ Xi are the second lightest neutralino and 
the lighter chargino). Also, we have not restricted ourselves to a supergravity (SUGRA) scenario [|l5[ 
but considered squark and slepton masses as free parameters within existing experimental bounds (this 
is a common practice in Tevatron SUSY analyses; see also Ref. [|lq]). 

The signals involving leptons get enhanced considerably if one remembers that the cascades leading 
to LSP pairs also produce single-or multi-leptons simultaneously. For example, heavier neutralinos 
(mainly X2) ^^d the lighter chargino (xi") can cause cascades where sleptons giving large contributions 
if they are so light as to be produced on-shell. On the whole, LSP pairs are in general produced in 
final states of the following types: 

X? X? + 0,1,2 or 3 ; + X with/ = e,/i,r (5) 

Charged leptons can be produced in the cascade, or in the decay Xi-pair; X denotes other states, 
possibly produced in association (jets, neutrinos). 

The sample result^ shown in Table 1 demonstrate that, of the different final states mentioned 
above, the 'pure jet + LSP pair' and the 'single lepton' final state have the highest rates in general, 

^Loop effects split the residual degeneration in the neutrino mass matrix; this permits to address the deficit of solar 
Ve M in the context of the present model (see also |P_0|). 

^Our results are cross-checked in the appropriate limits against, for example, the calculation in llq]. 
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Table 1: Lowest order (LO) cross sections in pb at y/s = 2 TeV for various final states and model parameters 
(/i' — —250 GeV). cr„'s (n — 0,1,2,3) are the cross sections for n-lepton + jets + Xi Xi + -9 (carried by 
neutrinos) final states; cr„; is the cross section for hadronically quiet n-lepton event. The numbers shown as .000 
are insignificant up to three decimal places. We used CTEQ-4L parametrization (from PDFLIB, |jl^) evaluated 
at Q = (mi+mj)/2 {rriij are the masses of the particles produced in the hard scattering). Next-to-LO corrections 
enhance ctq and ai by ~ 20% . 



the dominant channel being determined by whether a slepton is light enough to be produced on-shell. 
These rates are followed by that of the di-lepton final state. In the latter cases, however, our signals 
can receive significant contributions because, if one of the leptons (/i or r) pairs up with one of identical 
flavour and sign coming from LSP decay, the overall rate gets suppressed only by a single power of 
the LSP branching fraction rather than its square. Thus, by reconstructing one Xi and letting the 
other one decay in any allowed channel, one may get an enhancement by about a factor of 6 so long 
as it is enough to look for just one reconstructed W in the final state. The most convenient mode 
for reconstructing the W and hence, checking the mass-shell condition of Xi seems to be Xi — ^ ^^ 
(/ = /i,r) followed by W — > jj, since all the decay products are visible here. 

The standard model backgrounds to this kind of a final state {W + LSD) are quite suppressed. 
This is because the neutralino decay length in our case can be as large as between 0.1 mm and 1 
cm, [§, P when the decaying neutralino is of mass around 100 GeV or above. Such decay gaps result 
from the smallness of the R-parity violating coupling driving decays of the LSP. This coupling is 
determined by the basis-independent parameter vo which also controls the tree-level neutrino mass. 
The resultant displaced vertex- to which the W needs to be reconstructed- will distinguish our signals. 
For example, LSD backgrounds from tt-production followed by semileptonic decay of one of the bottom 
quarks produced in top-decay can be eliminated by proper identification of the primary and secondary 
vertices together with suitable isolation cuts on the leptons. Though most MSSM cascades potentially 
faking our signals are suppressed at some stage or other, there are processes like pp — > bb* followed 
by decays of the 6-squarks to a top, where both 1^'s and like-sign dileptons can be observed. However, 
the decay gap typical of our scenario can cause the latter to stand out. When R-parity is violated by 
trilinear terms, LSP's produced via the decay xt — ^ Xi^i followed by three-body leptonic decay of 
each LSP, can in principle give rise to LSD + W, where a similar decay gap as above can occur if 
the trilinear R-violating couplings are appropriately small. In order to avoid faking by signals of this 
kind, one could test whether one lepton in the LSD pair and the W originate from the same displaced 
vertex. In addition, here as well as for the ti backgrounds mentioned above, one can use the neutralino 
mass scale condition to suppress the backgrounds. 



Thus the finally suggested signals are of the form 

like-sign dimuons/ditaus + displaced vertex + a real W paired with one [iJt. (6) 

The condition of having one displaced vertex forces us to leave out those events where both leptons 
have their origins in cascades rather than in x? decays. 

In figures 1-3 we show some plots of the predicted number of signal events (of both like-sign dimuons 
and ditaus) expected at the Tevatron Run II (-^s = 2 TeV, with an integrated luminosity of 10 fb^^), 
once the above set of criteria are specified. The Higgsino mass parameter \x is taken to be —250 GeV 
in all the cases. No drastic effect in the nature of these curves are observed due to variations in [i. A 
common low-energy slepton mass has been assumed in each case. In addition, the physical stau masses 
are determined by left-right mixing depending on the value of tan/?. Each curve is truncated at the 
point where the lightest neutralino ceases to be lighter than all the sleptons. 

Let us discuss the most prominent features of these results: 
(i) For slepton masses of 150 GeV, leptonic final states can be observed via decays of on-shell sleptons 
produced in the preceding decays of neutralinos or charginos in the cascades. For low tan /3, like-sign 
dimuon and ditau events rates are closely comparable. As tan/3 increases, mixing in the stau sector 



lowers the mass of one physical state while still keeping it consistent with the current search limits [ 20 1. 
This results in larger event rates with ditaus and a consequent splitting between the dimuon and ditau 
curves as seen from figures 1, 2 and 3. Also, /xr-events with similar kinematic characteristics as our 
already described signals are expected in this scenario. The //r-type LSD event plots should correspond 
approximately to the sum of the dimuon and ditau curves for each combination of parameters, 
(ii) With the slepton masses on the higher side (225 GeV and 300 GeV), the event rates are mostly 
controlled by the all-jet channel for relatively low masses of the lightest neutralino. In these regions, for 
low tan /3, the numbers of muonic and tau events are very close together, and are relatively insensitive 
to the above variation in the slepton masses. However, with large enough gaugino masses, two-body 
decays of X2 ^-^^d Xi^ iii the cascades become possible, whereby rates for the single-and dilepton -|- 
x5 ~ pobir events get enhanced. This causes both the muonic and tau-events to go up comparably for 
lower values of tan /3 for which there is no substantial mass-splitting between the smuon and stau mass 
eigenstates. For larger values of tan/3 (tan/? = 10,30) in this region (see figures 2 and 3) there is a 
rather early onset of the two-body decays of X2 ^-'^d xi^ leading to a sharp enhancement in tauonic 
events while a relatively minor rise of this type is noticed for the corresponding muonic events at a 
larger value of gaugino masses. The difference is attributed primarily to the fact that for same tan /? 
the lighter r mass eigenstate is less massive than its muonic counterpart. This in turn implies that 
the two-body decays of gauginos involving muons open up for somewhat higher gaugino masses which 
suppress the signals from the production level as well. The rise in the events rates due to the above 
effects is offset by the fall in the rate of X2 ^^d Xi production when the masses of the latter (and 
therefore that of x?) go on increasing. That is why the curves showing the tau-rates as shooting up 
also show a fall following a peaking behaviour. 

(iii) However, for a slepton mass of 300 GeV these features are not clear from the graphs. The relevant 
cross sections surely get diminished due to heavier sleptons and LSPs in the processes leading to a 
smaller number of events. These, when plotted with the curves for lower slepton masses, fail to attain 
the required resolution due to obvious reasons. 

Our results indicate that rather copious numbers of uncut signal events^, unlikely to be faked by 
backgrounds, are expected over a large part of the parameter space. Comparing the upper and lower 
panels of each of of figures 1,2 and 3 it is seen that the number of events fall with a rise in the squark 
masses, since, as we have already mentioned, squark pairs form the main source of cascades here. Apart 
from that the general features of the events remain unaltered for fixed [i and tan /3. 

■^Actual analyses need to consider ID-efficiencies for /i, r and VK's, and isolation and pt cuts for /i, r. 



In fig. 4 we show the muon p^^-distributions in some sample cases. We have separately plotted 
the distributions coming from decays of the LSP {i.e. the ones emerging from displaced vertices) and 
those produced in earlier stages of the cascades (i.e. the ones that would be there even if R-parity 
were conserved). Increasing the LSP mass, the leptons from the cascades are less important for the 
overall signal strength, while the muons originating from the LSP tend to be harder. We have used a 
slepton mass of 150 GeV in these plots; the cascades, as demonstrated in Table 1, are considerably less 
significant for higher slepton masses. These conservative estimates show that a large number of muons 
in each case survive an p^-cut of the order of 15 GeV, usually adopted in Tevatron experiments [19j. 



Taking the overall detection efficiencies into account, we expect that well above 50 % of the events 
predicted in figures 1-3 can be salvaged after the experimental cuts employed at the Tevatron. 

In conclusion, the decays of the lightest neutralino can have quite distinctive signatures in the 
considered R-parity violating scenario. The characteristic signals are final states comprising of like-sign 
dimuons and ditaus, together with a real W, with a measurable decay gaps (0{lmm) or larger). These 
criteria eliminate any backgrounds, and distiguish the signal from those of other R-parity breaking 
models. An effective VF-reconstruction leads to observation of the lightest neutralino mass peak. 
These considerations motivate to a certain extent the urge for further refinement of the techniques of 
such reconstruction as well as of the resolution of displaced vertices at Run-II of the Tevatron. 

We thank V. Ravindran for useful discussions. 
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Figure 1: The predicted numbers of like-sign dimuon and ditau events (indicated by the labels /i and t on the 
respective curves) with a real W and a displaced vertex, as a function of mass of the lightest ncutralino (LSP) 
at the upgraded Tevatron with an integrated luminosity of 10 fb^^. Three values of the low energy slepton mass 
parameter (in GeV) have been used. Other parameters are : /i' ~ —250 GeV, tan/3 = 3. The degenerate squark 
mass rrig is 300 GeV in the upper panel, 400 GeV in the lower one. Notice the different vertical scale. 
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Figure 2: Same as in Figure 1, but with tan/3 = 10. 
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Figure 3: Same as in Figure 1, but with tan/3 = 30. 
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Figure 4: px distributions of muons for various values of tan/3. The upper panel shows the distributions for 
muons from a displaced vertex due to _R-parity violating decays of the lightest neutralino (the LSP). The bottom 
panel illustrates those for muons originating from usual MSSM cascades. The values of the MSSM parameters 

150 GeV, m^o = 100 GeV, fi' ^ -250 GeV. 
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